Characteristic frequency of the magnetic radiation of spinor condensates 
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The magnetic radiation of the fully-condensed states of ^^Na condensates have been studied. A 
narrow characteristic spectral line with a wave length cx "^^^ (A'' is the number of particles) and 
with a probability of transition oc N^''^^ emitted (absorbed) by the condensate was found. It 
implies that short wave radiation with a huge probability of transition can be obtained if numerous 
atoms are trapped. A new technique developed by the authors, namely, the analytical forms of the 
fractional parentage coefficients, was used to calculate analytically the matrix elements between the 
total spin-states. 
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The theoretical prediction and the subsequent experi- 
mental realization of scalar and spinor Bose-Einstein con- 
densation (BEC) are great achievements of 
scientific research. It implies that, under certain circum- 
stances, BEC as a new kind of aggregations of matter do 
exist. One can not rule out the possibility that they 
might exist in nature. If they do exist, an essential way 
to find them out is to measure the radiation emitted from 
them. For this purpose, one must know the character 
of the radiation. To the knowledge of the authors, this 
is a topic not yet been studied before. This paper is an 
attempt dedicated to the topic. 

For the scalar BEC of neutral atoms, neither charge 
nor magnetic moment are contained, therefore no radi- 
ation is involved. However, the spinor BEC contains 
magnetic moments, therefore magnetic multipole radia- 
tions may exist. For a single atom, the magnetic moment 
is small and the related radiation is weak. However, for 
the spinor BEC, the magnetic transition might be bene- 
fited by the Bose-enhancement, because numerous atoms 
are involved. Thus, an evaluation of the probability and 
wave length of the transition is meaningful. It would 
be interesting to see how large the particle number N is 
needed so that the magnetic radiation can be observed. 
This is an aim of this paper. For this purpose, the related 
initial and final eigenstates together with their eigenen- 
ergies must be known. 

A crucial point is the calculation of the matix elements 
for the transition between the total spin-states each con- 
tains numerous particles. When N is small, this can 
be achieved by using the fractional parentage coefficients 
(FPC) proposed about seventy years ago by Bachcr and 
Goudsmit.|6j] These coefficients are usually obtained from 
recursion relations. However, when N is large (say, 
N > 50), the procedure of recursion would become ex- 
tremely complicated. Therefore the FPC can not be used 
for many-body systems previously. Recently, we have 
succeeded to obtain the analytical forms of these coeffi- 



cients for spin-1 systems, thereby they can be used for 
marry-body systems as well disregarding how large N is 
[1, Thus, the main obstacle, namely, the difficulty 

in the calculation of related matrix elements, has been 
removed. 

Consider a condensate of N neutral spin-1 atoms 
trapped by parabolic and isotropic confinement so that 
both the orbital angular momentum L, total spin S, to- 
gether with their Z-components M and Ms are good 
quantum numbers. Let gij^^i and ^'^m s Ms 

the final and initial eigenstates, respectively. In gen- 
eral, the rank— J magnetic multipole transition proba- 
bility per unit time reads [Tl] 
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where uj is the frequency of the emitted (absorbed) pho- 
ton, fiatom is the magnetic moment of the atom. 
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where j denotes the j — th atom, Fj the spin operator, 
YjK the spherical harmonic, K ranges from — J to J. 

Let the components of be Fjx, Fjy, and Fj^, and 
Fj,±i = -^{Fjx ± iFjy) , F.jfl = Fj^. For the dipole 

transition (Ml), Si,±i,j = ipy^F,-^ ^i, Si,o,j = 

|:Fj.o- For the quadruple transition (M2), 'E.2,±2. j — 
±\/TOrjyi^zpiFj, zpi, S2,±i,j = T\/5rj(YioFj^ ^1 ± 

yi,Ti-fj\o), ^2,0, J = -y/§^j(^ii-Pj" -1 -yi,-iFj,+i + 
2Yi,i)Fjfi). For J > 3, the transition probabilities would 
be much smaller and therefore omitted. 
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Obviously, the Ml operator keeps all the good quan- 
tum numbers completely unchanged, except Mg. If 
spin-orbital coupling is not considered, and if there is 
no external magnetic field, the Ml-transition can be ne- 
glected. Therefore, we shall focus on the M2-transition. 
During the transition, \L' - L\ < I, \S' - S\ < I and 
M' + M'g + K = M + Ms are required. In particular, 
the parity must be changed, 

Let TO be the mass of the boson and Wo be the fre- 
quency of the parabolic confinement. When Hwo and 
yjh/mwo are used as units, the Hamiltonian reads 



(5) 
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where Ui-j = 5{ri~rj)Oij, 0,j = (co -I- C2Fi • F^). 
this paper we consider only the transitions related to the 
fully condensed (FC) states 



^SMs 



SMs 



(6) 



where '&^sm^ is an all-symmetric normalized total spin- 
states with the good quantum numbers S and Ms- S is 
allowed to range from N, N — 2, iV — 4, • • •, to 1 or O.fl^. 
[Tsl The single particle state (j)s{r) = us{r)/\/4^r has 
orbital angular momentum zero, and it depends on S in 
general. The spatial wavcfunction (t>s{r) can be obtained 
by solving the generalized Gross-Pitaevskii equation 
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and £5 is the chemical potential. The eigenenergy of a 
FC-state depends surely on S, it reads 



where N — 1 particles condense into the same (figir) = 
-Ug{r) , while a particle is excited to a higher state 

^Sa'/s ^ normalized total 
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^M,s(^) = T4(0>^iAf(r) 

spin-state with the good quantum numbers S and Ms, 
and with A = |7V— 1, 1}, where S is allowed to range from 
1 to iV — 1. [all^l Obviously, the summation over i in 
(10) assures that "^m^sms all-symmetric as required. 



and Us satisfy the set of coupled equations 
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where h, ^ + ^ + r\ = {&^sili's\0 ,i^\Q^s^s) 

{^^J^k), and = {efjIillO^A^f^il + ©sm^)- By 



using the FPC, gg and gg can be derived as| 
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Since g^ and g| are known, eigenstates can be obtained 
by solving (11) and (12), the associated eigenenergy reads 
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(9) where hi^ J dr u%hiu% , hl^ J dr u^shou^ , Tf 



where ho= J dr ushous and T5 = ^ ^Wsl"^- Since 
Eg depends on 5, all the FC-states with S ranged from 
N, to 1 or spread into a band, the FC-band.[12] 

It is obvious that the M2-operators can excite only 
one p— wave particle. In other words, a FC-state would 
transit to a final state with one particle excited. From the 
rule of outer product {l}(g>{iV- 1} = {N} + {N - 1,1}, 
the spatial wave functions of the final states may have two 
choices of permutation symmetry A. The choice A = {N} 
is just an excitation of the center of mass (cm.), the 
related photon energy is just hwc,. This is a trivial 
case with a very low photon energy. For the choice 
A — {N — 1, 1}, the final state with an excited particle 
can be written as M 
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From now on let S (5') denotes the total spin of the 
initial (final) state. For M2 transitions we have S' = 
S or 5 ± 1. From (15) and (9), the energy difference 
{AE)s's = Efcjf ~ E% can be known ( \[l\E)s's\ is the 
photon energy). 

By solving eq.(7), (11), and (12), not only the photon 
energies can be known, the probability of M2-transition 
can also be obtained. To compare with experimental 
observation, we define 
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By inserting (6) and (10) into (3), we have 
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where Is's = J dr u^g,rus and 

Ms 

+i(0L'^,'4-ii^i,-ii'?L'Il>P} (18) 

It turns out that the quantity inside the brackets in 
(18) does not depend on Ms- By using the FPC and 
refer to the appendix of [f^, we have 

Ys's^{'iS+l){N~S){N + S+l)/N^ (19) 

Ys,s^{N-S){S + 2)lN (20) 

and 

Ys's^iN + S+l){S-l)/N (21) 

if S" = 5, 5 + 1, andS — 1, respectively. 

The numerical results of Na atoms are presented in 
the follows as examples. The interaction of Na has 
Co = 6.774 X 10-*v^ and C2 = 2.117 x IQ-""^, 
Wo — 300 sec" is given in this paper. {AE)s's are given 
in Fig.l. There are three points noticeable, (i) For the 
S to S' = S transition, there is a turning point appear- 
ing at S/N = l/-\/3 . This is associated with a sudden 
change of the structure of the excited particle as shown 
later, (ii) {AE) 5/ 5 can be either positive or negative. It 
implies that, when N is very large, the FC-states may 
not be the lowest, (iii) The magnitude of {AE)s's of the 
5 to 5" = S* — 1 transition is relatively small, thus the 
probability of this transition is negligible. 

The probability Pg,'g is given in Fig. 2, where only the 
curves for S' — S transition (solid line) and for 5" — S+1 
transition (dashed line) can be seen. There is a striking 
sharp peak at S/N — 1/^/3 for the S' — S transition 
(just at the turning point of {AE)s's ), the height of 
the peak is 1.44 x 10~'^'*/sec, it is so high that only the 
bottom of the peak was seen in the figure. Accordingly, 
there would be a narrow spectral line emerged more than 
ten times brighter than the one contributed by the lower 
peak. The lower peak is much lower and much broader, 
it is peaked at S/N = 0.96 and is contributed by the S' — 
S+1 transitions. From Fig.l we know that the narrow 
peak has photon energy 127.9, while the broad peak has 
photon energies ranging from about 220 to 380. Thus, the 
feature of the spectrum is clear, namely, a bright narrow 
line together with a background at its violet side. The 
narrow line is called the characteristic line of the system 
thereafter. 

When N varies, the features of Fig.l and 2 remain un- 
changed, but the magnitudes of {AE)s's and Ps's are 
changed remarkably. In Fig. 3 the wave lengths Xs's as- 
sociated with the photon energies of the narrow peak and 
the lower peak against log{N) are plotted. Where the two 
parallel straight lines imply that Xs'S are proportional to 



N in both cases, thus very short wave radiation can 
be obtained if N is sufficiently large. In Fig. 4 the proba- 
bilities Pgfg are plotted. The solid line as a straight line 
implies that the probability of transition of the charac- 
teristic line is proportional to iV^^/^. The power 17/5 
is striking, it implies that the probability can be terribly 
large when N is very large. The dashed line implies that 
the probability associated with the S' = S +1 transition 
is much lower, and it increases with A'' much slower. For 
a quantitative example, if A > 10^'*, Psfg of the charac- 
teristic line would be > lO^'sec^^ and Xs'S < 150cm. 

Let us study the physical background of the above find- 
ings. When N is very large, we can neglected the second 
term of (11), then Ug, can be obtained by using further 
the Thomas- Fermi approximation, [l^l (neglecting the ki- 
netic energy). When the Ug, obtained in this way is in- 
serted into (12), the equation can be rewritten as 

ll(-^ + ^) + UefM^^s,4, (22) 

Where 

C^e// = ^(r^-r2)C(r,-r)-K^ (23) 

where — {15Ngg, /Aw)^^'^ and C(?'a — r) = 1 if r < Ta , 
or = otherwise. 

It turns out that the feature of Uef / is extremely sen- 
sitive to gg, / gg, . When — S" is even and S' is close to 
N/V5 (or N-S' is odd and S' is close to A^), fflVffs, « 1, 
thereby Ueff is nearly a constant in the broad domain 
r < Tq, but has a sharp border at as shown by the 
S' = /-s/S line of Fig. 5. In this special case, Ug, spreads 
widely from to Tq, and results in having a very large 
Is'S, and thereby a very large probability of transition. 
Alternatively, when ggi/ggi < l,C/e// is attractive as 
shown by the S' = N curve of Fig. 5. In this case the 
excited particle would be distributed close to the cen- 
ter. When g^i/ggi > 1, Ueff is repulsive as shown by the 
S" = curve of Fig. 5. In this case the excited particle 
would be distributed very close to the surface of the con- 
densate. To give quantitative data related to Fig. 5, Ug/ is 
mainly distributed in the domain < r < 11 if S" = A^, 
< r < 148.5 if S' = A^/V3, and 145.5 <r < 148.5 if 
S' = 0. These data exhibit the great difference in the 
wave functions of the excited particle, which results in 
the dramatic variation of Is's- When S' = S we have 
Is's = 4.24, 14.07, 87.21, 0.81, and 0.04 if S' = 0, 
57Af/100, N/VS, 58A^/100, and A^, respectively. These 
data exhibit that, when S' varies across N/\/3, a great 
increase and a great decrease of Is's occur successively, 
this is the origin of the narrow peak shown in Fig. 2. 

For the transition with N-S' odd and S' = S+1, 
a similar increase of Is's occurs if S" is approaching A^. 
However, in this case Ys's tends to zero (refer to eq.(20)). 
As a result of the competition of Is's and Ys's, the lower 
and broad peak appears adjacent to S" = A^ as shown in 
Fig.2. 
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In summary, wc have found that the spin-dependent 
interaction is crucial to the magnetic radiation of spinor 
BEC. The effect of the interaction is embodied by the ra- 
tio Qsi/gs', which affects Ugff seriously. Accordingly, the 
excited particle has three choices of structure, namely, 
being close to the center, close to the outer surface, and 
widely distributed from the center to the surface. The 
wide distribution leads to a great increase of -P/Zg. How- 
ever, this case occurs only if S' « N/\/3 and S' = S, or 
S' N and S' = S + 1. This leads to the appearance of 
the characteristic line together with a broad background 



as shown in Fig. 2. The N-dependence of Xs's and Pgtg 
has been studied, and a N~'^/^ and N^"^/^ dependence, 
respectively, was found. This implies that short-wave 
radiation with a large probability of transition can be 
achieved if numerous particles are trapped. 

It is mentioned that only the magnetic radiation of the 
FC-states has been studied in this paper, it is still a long 
way to understand the spectroscopy of the spinor BEC. 
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FIG. 1: {AE)s's = Ef' - E's against S/N for ^^Na atoms 
with AT = 5 X 10^^ The solid curve is for S' — S . the dashed 
curve is for S' — S+1, the dash-dot-dot curve is for S' = S—1. 
Wo = 300 sec^^, the unit of energy is hwo- The notations and 
the value of Wo are the same in the following figures. 




S/N 



FIG. 2: The probability P/;° of the M2-transition of the FC- 
states of ^^Na against S/N. iV = 5 x 10^^ . 
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FIG. 3: The wave lengths Xs's (in cm) of the radiation of the 
FC-states against log{N). The solid line is for the character- 
istic line which has 5* = and S' — S, the dashed line is 
for the lower peak which has 5" = 96Af/100 and 5" = 5" + 1. 
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FIG. 4: The probability of transition (in sec ^) against 
log{N). Refer to Fig.3 . 




FIG. 5: Ueff{r) for the excited particle with three cases of 
S' , N — S' is assume d to be even. The units of energy and 
length are hwo and \fhfrrvw'o., respectively. 



